Collective modes of trapped Fermi gases with in-medium interaction 
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Due to Pauli blocking of intermediate states, the scattering matrix (or T matrix) of two fermionic 
atoms in a Fermi gas becomes different from that of two atoms in free space. This effect becomes 
particularly important near a Feshbach resonance, where the interaction in free space is very strong 
but becomes effectively suppressed in the medium. We calculate the in-medium T matrix in lad- 
der approximation and study its effects on the properties of collective modes of a trapped gas in 
the normal-fluid phase. We introduce the in-medium interaction on both sides of the Boltzmann 
equation, namely in the calculation of the mean field and in the calculation of the collision rate. 
This allows us to explain the observed upward shift of the frequency of the quadrupole mode in the 
collisionless regime. By including the mean field, we also improve considerably the agreement with 
the measured temperature dependence of frequency and damping rate of the scissors mode, whereas 
the use of the in-medium cross section deteriorates the description, in agreement with previous work. 

PACS numbers: 67.85.Lm 



I. INTRODUCTION 

In the last few years, experiments on collective modes 
in ultracold trapped Fermigases at Duke University [l], Q 
and at Innsbruck [H, 0, [f| [y, 0] revealed a lot of interest- 
ing information about the different regimes which can be 
realized in these systems. This is mainly due to the pos- 
sibility to change not only the temperature of the gas, 
but also the interaction between the atoms by adjust- 
ing the magnetic field in the vicinity of a Feshbach res- 
onance. One of the objectives of these experiments was 
to detect the transition from the normal to the super- 
fluid phase. However, although the frequencies of the 
collective modes arc generally different in the hydrody- 
namic and in the collisionless regime, it was recognized 
that a hydrodynamic behaviour of the gas is not an un- 
ambiguous sign of superfluidity. In fact, hydrodynamic 
behavior can be a consequence of superfluidity, but also 
of a sufficiently high collision rate in the normal-fluid 
phase [8j. Recently it was therefore proposed to distin- 
guish the superfluid, the collisionally hydrodynamic, and 
the collisionless regime @. 

The experimental results suggest that the transition 
from the superfluid to any of the normal-fluid regimes 
is always accompanied by a very strong damping of the 
collective modes. While the frequencies of the collec- 
tive modes in the zero-temperature limit can easily be 
predicted in the framework of superfluid hydrodynam- 
ics, literature on the finite-temperature case is relatively 
sparse. One suggestion put forward by Griffin et al. was 
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to apply Landau's two-fluid hydrodynamics @, but this 
requires that the collision rate in the normal-fluid compo- 
nent is sufficiently high, i.e., much higher than the mode 
frequency. This is usually not the case, since in a trapped 
system the collective modes have frequencies which are 
at least of the order of the trap frequency. In the op- 
posite limit, i.e., if the normal-fluid component is colli- 
sionless, as it should be in the weakly interacting limit, 
a quasiparticle transport theory which couples superfluid 
hydrodynamics to a Vlasov equation for the normal com- 
ponent was suggested by one of the authors [Ty 1 , EH, E2] • 
Although this theory provides a qualitative explanation 
of the damping of the collective modes near the transi- 
tion from the superfluid to the collisionless normal-fluid 
phase, it cannot give a quantitative description of the 
recent experiments since these are generally done in the 
strongly interacting regime (near a Feshbach resonance). 

In the present paper, we will consider the case of at- 
tractive interaction, i.e., negative scattering length a < 0, 
from the weakly interacting regime far from resonance up 
to the unitary limit, |o| — > oo. We will concentrate on 
the normal- fluid phase, i.e., on temperatures above the 
critical temperature Tq of the superfluid-normal phase 
transition. The experimental data taken at these tem- 
peratures still show a rather strong damping of the col- 
lective modes, since the experiments are in fact neither 
in the collisionally hydrodynamic nor in the collisionless 
limit. This intermediate regime has been studied theoret- 
ically in the framework of the Boltzmann equation (some- 
times also called Boltzmann- Vlasov or Landau- Vlasov 
equation). In the framework of the Boltzmann equa- 
tion, interactions between the atoms enter in two dif- 
ferent places: On the one hand, each particle moves in a 
mean field characterizing the average of the interaction 
with all other particles. On the other hand, the parti- 
cles undergo two-body collisions, which are determined 
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by the scattering cross-section. 

In some of the early literature on this subject [H, EE 
EE EE E3 ' tne mean field effects were considered within 
the Hartree approximation. However, the most recent 
work by Bruun et al. @, EH concentrates mainly on the 
strongly interacting regime, where the Hartree approxi- 
mation breaks down: While more and more sophisticated 
models of the collision term were developed, mean field 
effects were completely neglected @, EH- Therefore, the 
predicted frequencies in the collisionless limit are those 
of an ideal gas, and in the hydrodynamic limit, the equa- 
tion of state is that of an ideal gas, too. However, at least 
the experiment of Ref. [4 : ] clearly shows the importance 
of the mean field shift in the collisionless regime. 

The aim of the present paper is to include mean field 
like effects on the propagation of the particles into the 
Boltzmann equation in a way which is appropriate also 
for the strongly interacting case. Our approach is based 
on the in-medium scattering amplitude (T matrix), cal- 
culated in ladder approximation. This allows for a unified 
description of mean field like effects on the propagation of 
the particles (self-energy) and the modified cross-section 
entering the collision term of the Boltzmann equation. 
We point out that in order to have a consistent theory, 
it is important to use in the description of the collective 
modes, which are small variations around equilibrium, 
the same self-energy as in the calculation of the equilib- 
rium density profile. This can be seen, e.g., by looking 
at the frequency of the sloshing mode, whose frequency, 
according to the Kohn theorem S] , must be equal 
to the corresponding trap frequency. 

Our paper is organized as follows. In Sec. UH we con- 
centrate on the description of the system in equilibrium. 
After illustrating the breakdown of the Hartree approxi- 
mation, we review the treatment of the interacting Fermi 
gas in the framework of the ladder approximation, intro- 
ducing the in-medium scattering amplitude, the single- 
particle self-energy, and the in-medium cross-section. We 
discuss the quasiparticle approximation which is neces- 
sary for describing the system in the framework of the 
Boltzmann equation and which allows us to define a 
quantity similar to the mean field in Hartree approxima- 
tion. Wc show the resulting density profiles and discuss 
the limits of validity of the quasiparticle approximation. 
In Sec. IIII1 we turn to the description of collective modes 
in the framework of the Boltzmann equation, which now 
includes the "mean field" and the in-medium cross sec- 
tion. Using the standard method of taking moments of 
the Boltzmann equation in phase space, we obtain semi- 
analytic expressions for the collective mode frequencies 
which depend only on integrals of equilibrium quanti- 
ties. In Sec. IIV1 we discuss our numerical results for 
the scissors, radial quadrupole, and breathing modes and 
compare them with the available experimental data for 
temperatures above Tq- Finally, in Sec. [Vj we draw our 
conclusions and give an outlook to further developments. 

Throughout the paper, we will use in the derivations 
units with h — ks = 1 (H = reduced Planck constant, fee 



Boltzmann constant). 



II. IN-MEDIUM EFFECTIVE INTERACTION 

A. Breakdown of the Hartree approximation 

Let us consider a two-component (f, J,) uniform Fermi 
gas in the normal phase. As long as the range of the inter- 
action is small compared with the mean distance between 
atoms, one can assume a zero-range (5 function) interac- 
tion between atoms of opposite spin, and the hamiltonian 
reads (in second quantization) 



H = 



d 3 i 



V 2 \ 

, (i) 

2m - 1 1 / 



where m and g denote the atom mass and the coupling 
constant and ip is the fermion field operator. We as- 
sume that the interaction is attractive, i.e., g < 0. The 
coupling constant is related to the atom-atom scattering 
length a via 



fJ 
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In Hartree approximation, the single-particle energies of 
t and I particles are shifted by Uhi — gp\ and Uhi — 
gp\i respectively. The exchange or Fock term vanishes 
since the interaction is only between atoms with opposite 
spin. From now on, we will assume that both spin states 
are equally populated, and we denote by p = p-\ = Pi 
the density per spin state. Then the Hartree shift is the 
same for both spin states, and we may write 



U H = gp- 



(3) 



Let us now calculate the density as a function of the 
chemical potential p. For our purposes it is enough to 
consider the zero-temperature limit, i.e., temperatures 
T <C ep i where 



k 2 

ep = - — and 

2 m 



kp = (67r 2 p) 



1/3 



(4) 



are the Fermi energy and Fermi momentum, respectively. 
[Note that Eq. ^ defines ep and kp for a uniform sys- 
tem, i.e., it remains valid even if the occupation numbers 
do not resemble a step function because of temperature 
or correlation effects. For a trapped system, however, we 
will use a different definition of ep and kp, see Sec. Ill F[ ] 
At zero temperature, the relation between ep and p is 
given by eF = p — Uh ■ Substituting Eq. ([3]) into this 
relation, one obtains the following cubic equation for kp: 
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It is easy to see that this equation does not have a so- 
lution if p exceeds a critical value given by p ma x = 
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7r 2 /(24ma 2 ), corresponding to a maximum density of 
p m ,n, T . = 7r /(48|a| 3 ). The same value was found in Ref. 
[211 ] as the density where the system becomes unstable 
against separation into a low-density (gas) and a high- 
density (solid) phase. 

If the above arguments were correct, a low- 
temperature Fermi gas with attractive interaction should 
be unstable as soon as kp\a\ > ir/2. However, we know 
from experiments that ultracold Fermi gases are stable 
throughout the BCS-BEC crossover, including the uni- 
tary limit kp\a\ — > 00, because the system prefers to form 
pairs instead of separating into two phases [22I ] . The in- 
stability is simply an artefact of the Hartree approxima- 
tion and not physical. 



B. 



T matrix 



In the Hartree approximation as described in the pre- 
vious subsection, the coupling constant g was related to 
the scattering length a in free space. This means that it 
was implicitly assumed that the scattering amplitude it- 
self is the same in the gas as in free space. As we will see, 
this assumption is the origin of the unphysical instability 
of the Hartree approximation at high density or strong 
interaction. For instance, as pointed out in Ref. [23| the 
scattering amplitude becomes proportional to 1/kp in- 
stead of a at high density. 

The approximation scheme we adopt here in order to 
calculate the in-medium scattering amplitude is based on 
the non self-consistent T matrix approximation. In this 
approximation, the T matrix is given by the resumma- 
tion of ladder diagrams, and it depends only on the total 
energy E = cu + 2/i and the total momentum k of the two 
atoms: 



r(w,k) = 



l- s j(w,k) 



(6) 



where J denotes the non-interacting two-particle Green's 
function. Within the imaginary-time (Matsubara) for- 
malism [24j, the latter is given by 

J(iu; N ,k) = -T J -±jL 0oK,k/2-q) 

n odd 

x g Q (uj N - uj n ,k/2 + q) , (7) 

where ujn and w n are, respectively, bosonic and fermionic 
Matsubara frequencies, and Qo{w n ,\s) = l/(ioj n — ££) is 
the free (Matsubara) Green's function, — k 2 /(2m) — fi 
being the free single-particle energy. After evaluation of 
the sum over n, the retarded function J(w, k) is obtained 
as usual by analytic continuation. The result reads 



J(w, k) 



d 3 q 



1 n k/2+q n k/2-q 



k/2+q Mc/2-q 



IT) 



(8) 



where n° = l/[exp(/?£°) + 1], with (3 = 1/T. 



The problem with Eqs. ([I])-© is that J is divergent. 
To resolve this problem, one can introduce a momentum 
cut-off A, determine the coupling constant g as a function 
of A such that one recovers the correct scattering length 
a in free space, and finally take the limit A — > 00 keeping 
the free-space scattering length a constant [25[. In this 
way one obtains for the T matrix in free space 



r (£,k) 
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m 1 + iaq c 
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where q cm — ^JmE — fc 2 /4 is the on-shell momentum in 
the center-of-mass (CM) frame. If we now decompose J 
into the two-particle Green's function in free space, J , 
and a medium correction, J, such that J = Jq + J , we 
may write 



r( w ,k) 



Ana 
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m 1 + iaq c 



A-na j 
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with q cm = 'm(uj + 2fi) — fc 2 /4. Even without cut-off, 
the medium contribution J is finite and given by 



J(w,k) = 



d 3 q 



7° 

l k/2+q 



1° 

'k/2-q 



(2vT) 3 W 



^k/2H 



?k/2-c 
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The imaginary part of J can be given in closed form: 



Im J(u>, k) 



m T . 
2-ixk 



■In 



(12) 



where ^ = (k/2±q cm ) 2 /(2m) — fi. The real part is then 
computed numerically via a dispersion relation, 



Re J(w, k) = -V 



du' Im J(w',k) 



(13) 



As a by-product, the in-medium T matrix allows us 
to determine the critical temperature Tq of the system, 
i.e., the temperature below which the system becomes 
superfluid. As realized by Nozieres and Schmitt-Rink in 
their pioneering paper [22| , the Thouless criterion which 
relates Tc to the temperature where the T matrix devel- 
ops a pole at the Fermi level (i.e., at w = 0), remains 
true at all couplings. Since the pole always appears first 
at k = 0, the critical temperature can be obtained from 
the equation 



Rc J(u) = 0, k = 0; T = T c ) = - — . 

Ana 



C. Self-energy 



(14) 



Contrary to the zero-range interaction used in Sec. 
Ill A[ the in-medium vertex function T is now momentum 
and energy dependent. This complicates the calculation 
of the single-particle energy shift. Even the concept of 
such an energy shift may be questioned if there are no 
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well denned quasiparticles, as it is the case in the "pseu- 
dogap regime" [2y|. In any case, the appropriate object 
to calculate is the single-particle self-energy E, which is 
well-defined and does not rely on the existence of quasi- 
particles. The ladder self-energy can be written within 
the Matsubara formalism as 

d 3 p 



S(iw n ,k)= ( -^r^T Go(uJ„>,k) 

J [Z7T > n> odd 



x T(ioj n + iuv,p + k) . (15) 

Using analytic continuation to real energies, we find for 
the imaginary part of the retarded self-energy: 



Im E(w, k) 



d 3 p 



1 



- 1 



xlmr( w +^k + p), (16) 

which has to be evaluated numerically and from which 
the real part can be obtained by a dispersion relation 
analogous to Eq. (fT3"|) . 

Since we are going to use the Boltzmann equation for 
the description of collective modes, we are implicitly as- 
suming that the quasi-particles, especially near the Fermi 
surface, are well defined, which is of course a limitation 
of the range of applicability of our approach [2^] . If there 
are well-defined quasiparticles, this means that their dis- 
persion relation £k, determined by the poles of the single- 
particle Green's function, can be obtained from 



Re£(£ k ,k) 



(17) 



Such a treatment is probably desirable but beyond the 
scope of the present work. Here, we will completely 
neglect any energy and momentum dependence of the 
self-energy. Since we are mainly interested in momenta 
around the Fermi momentum, and hence energies around 
the Fermi energy, we will make the approximation 

£ k -£k + ^ ! with U = ReE(0,k fl ) ) (18) 



where k^ = ^/2mp [which can actually be quite differ- 
ent from kp as defined in Eq. Q]. It is the quantity U 
which will take the role of the mean field potential in the 
Boltzmann equation. 



D. Density 

The fact that the particles are interacting among each 
other changes strongly the equation of state of the sys- 
tem, i.e., the relation between the chemical potential p 
and the density p. Within the Hartree approximation, 
the relation p(p) can trivially be obtained, whereas in 
ladder approximation, the calculation of the density for 
a given chemical potential p is more involved. 

In principle, the density can be obtained as 



d 3 p 
(27p 



n odd 



T V £K,p) 



(19) 



where Q denotes the full single-particle Green's function, 
which according to Dyson's equation is given by Q^ 1 = 
Qq 1 — £ . Here we will restrict ourselves to an expansion 
of Q up to first order in S, i.e., we write 

g ~ g + GoEGo . (20) 

This leads us to the following expression for the density: 
P = Po+Pi, (21) 
where the free (uncorrelated) part, po, is given by 
d 3 P _ \ - „ f d 3 p 



Po 



J £fr T X Igfr*- < 22 ) 

n odd 



while, after a lengthy calculation, the expression for the 
correlated part, p\, can be written as 

Pi = J J^y T £ ^o( w n,k)E(iw„,k) 



odd 

d 3 k f du 1 



(2tt) 3 J 2tt -I dp 



Imln[-r(a;,k)]. (23) 



This is actually the result for the density initially given 
by Nozieres and Schmitt-Rink (NSR) in Ref. ^ and by 
Sa de Melo, Randeria, and Engelbrecht (29[. 

Unfortunately, the density formula given above is not 
suitable for being used as the ground-state density in the 
Boltzmann equation. In the Boltzmann equation, the 
density is expressed as an integral over the distribution 
function /: 



d 3 P f(p) 



(24) 



and / must reduce to a simple Fermi function. Our quasi- 
particle(QP) approximation for the density consists in 
using the single-particle energy £ p as defined by Eq. (fT5"|) 
in the calculation of the occupation numbers, i.e., 



feq(p) 



XI A 



1 



(25) 



where we have introduced the abbreviation A = (27r) 3 . 
In this way we obtain an alternative method for calcu- 
lating the density as a function of p. Fortunately, it 
turns out that in most cases both ways of calculating 
the density give similar results. In order to demonstrate 
this, we show in Fig. [T]the dimensionless ratio p/ep as 
a function of the parameter l/(kpa) defining the inter- 
action strength for various temperatures (in units of the 
Fermi energy). One can see that the interactions lead 
to a reduction of p with respect to the ideal gas result. 
At weak interactions, NSR and QP curves are in perfect 
agreement, but in the strongly interacting case there can 
be noticeable differences, especially at low temperatures 
close to Tc, indicating the breakdown of the QP approx- 
imation. This is not surprising, since in this case the 
system is in the "pseudogap regime" [26j . 
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FIG. 1: The ratio n/ep, with /i calcultated with the NSR 
density and the QP density, is plotted as a function of — l/fc^a 
for different temperatures: T/Tf = 0.05, T/T F = 0.3 and 
T/T F = 0.5. The curve for T/T F = 0.05 stops at l/k F a » 
1.36, since at this interaction strength the critical temperature 
is reached. 



E. Cross section 

The interaction between atoms is not only responsible 
for the single-particle energy shift. It also determines the 
collision rate of the atoms, which will play a central role 
for the properties of collective modes. The important 
quantity is the cross section a. In the case of a zero- 
range s-wave interaction, the cross section for two atoms 
in free space with momenta p 1 and p 2 before the collision 
and p'j and p' 2 after the collision is given by *2T] 



1 + {qaf 



(26) 



where q = (p 1 — p 2 )/2 is the incoming momentum in the 
CM frame and f2 is the solid angle after the scattering, 
i.e., dU. = 2iTsm9d9, where 9 is the angle between q and 
q' = (p[ — p 2 )/2 (note that |q| = |q'| because of energy 
and momentum conservation). 

As pointed out in Refs. @, [28|, the cross section is 
strongly modified by medium effects. In terms of the T 
matrix, the in-medium cross section can be written as 



da 

dn 



m / k 2 
47T V4m 



m 



-2/x,k 



(27) 



which now depends also on the total momentum k = 
Pi + P2 = Pi + P2 °f the two atoms. 

It is not easy to see a priori what will be the effect of 
these medium modifications, since depending on k and 
q the cross section can be enhanced or reduced as com- 
pared with the cross section in free space. This is shown 
in Fig. [2] where we display the cross section as a function 
of q for the case k — (where the medium effect is sup- 
posed to be strongest) for various temperatures for fixed 
interaction strength kpa = — 1. The strong enhancement 
of the cross-section near the critical temperature is a pre- 
cursor of the singularity of the T matrix at the critical 
temperature @, 
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FIG. 2: The ratio of the in-medium cross section and the vac- 
uum one is displayed for total momentum k = as a function 
of the relative momentum q. The results are shown for differ- 
ent temperatures for fixed interaction strength fen = — 1. 



Let us mention that the same effect has already been 
found some years ago in the in-medium nucleon-nuclcon 
cross section in low-temperature nuclear matter [3fj| . 



F. Local-density approximation 

Until now we considered a uniform system where the 
atoms are not trapped by an external potential. In or- 
der to include the trap potential Vr(v), we will make 
use of the local-density approximation (LDA), where the 
system is treated as locally homogeneous, with an r de- 
pendent chemical potential which is given by 



Mo 



V T {r). 



(28) 



This approximation should be valid as long as the po- 
tential varies only slowly, i.e., on length scales which are 
large compared with 1/kw- By the way, this condition is 
also necessary for the validity of the Boltzmann equation 
which will be used later to describe the collective modes. 

Within the local-density approximation, all equilib- 
rium quantities of the system discussed in the preceding 
sections, like T, U, p, da/dft, etc., acquire an additional 
r dependence via the dependence of /i on r. 

In practical calculations, we will use a harmonic trap 
potential 



V T {r) 



rn 
~2 



E < 

t—x,y,z 



(29) 



Besides the fact that experimental traps are almost har- 
monic, the use of a harmonic potential has the advantage 
that it is sufficient to calculate the equilibrium quantities 
once for a spherical trap with the same number of atoms 
and the average frequency lu — (u> x tjj y u) z ) 1 / 3 . Then, the 
equilibrium quantities in the deformed trap can easily 
be obtained from the corresponding ones in the spher- 
ical trap by the change of variables fj = riUi/uj (for 
i = x,y,z). 
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FIG. 3: Density profiles for 400000 atoms with l/k F a = 
—0.45 at various temperatures, calculated using the NSR for- 
mula (solid lines) or the QP approximation (long dashes). 
For comparison, the corresponding density profiles of an ideal 
Fermi gas (short dashes) are shown, too. The length unit is 
the harmonic-oscillator length lho = \J\Jrrvjj. 

As an example, we show in Fig. [3] the density profile 
of N = 400000 atoms in a trap with fixed interaction 
strength 1/kpa = —0.45 at different temperatures. Note 
that we follow here the usual convention of the experi- 
mental papers, using for the trapped gas a definition of 
kp and ep which is different from that used in the uni- 
form case [cf. Eq. Q]. In the trapped case, Uf and cf 
refer to the values of fcj? and in the center of the trap 
calculated for an ideal Fermi gas at zero temperature. 
Therefore, ep and hp are determined by the number of 
atoms, iV, and the average trap frequency as follows: 

e F = (3N) 1/3 u>, k F ^V2r7^. (30) 

One can see that the interaction leads to a noticeable 
change of the density profile, especially at low temper- 
ature. The agreement between the NSR and QP re- 
sults is almost perfect in the case T/Tp — 0.5, while at 
T/Tp = 0.21, corresponding to the critical temperature, 
the QP density is too small, indicating again the break- 
down of the QP approximation in the pseudogap regime. 
However, it should be mentioned that for 1/kpa = —0.45, 
one has to go very close to Tc in order to see this effect; 
it is more important in the unitary limit. Let us mention 
that similar density profiles can be found in the litera- 
ture, even for temperatures below the critical one [3l| . 

III. COLLECTIVE MODES 

A. Linearized Boltzmann equation 

After the discussion of static properties, let us now 
turn to the description of collective modes of a trapped 
Fermi gas. We remind the reader of three assumptions 
mentioned earlier: (a) the density is supposed to vary 
(due to the trap as well as the collective motion) only on 



large length scales; (b) the temperature has to be above 
the superfluid transition temperature Tc; (c) the quasi- 
particles near the Fermi surface have to be well defined. 
In addition, the Boltzmann equation is only valid if (d) 
the time dependence of the excitations under consider- 
ation (in the case of the collective modes under consid- 
eration this time scale is set by the the trap frequency) 
are slow compared to the "correlation time" [32J. Un- 
der these assumptions, the dynamics of the system can 
be described by the semi-classical distribution function 
/(r, p, t) whose time evolution is governed by the Boltz- 
mann equation (33j: 

/ + r • V,,/ + p • V p f = -/[/] , (31) 

where r and p satisfy the classical equations of motion. 
In the case of a quasiparticle dispersion relation as given 
by Eq. (jTSJ) with an r dependent chemical potential as 
given by Eq. ([28]) , the velocity and acceleration read 

r = V p £p = ^ , (32) 
m 

P = -V r £ P = -V r (V T + U) , (33) 

since, within LDA, £ p -> £ p (r) = p 2 /(2m) + V T (r) + 
U(r) — no. Note that there are two sources of r depen- 
dence of the self-energy U. In equilibrium, U depends on 
r only via /i(r) = fj,Q — Vr(r). More generally, in partic- 
ular out of equilibrium, the self-energy depends on the 
distribution function /, i.e., we may write U = U[f]. 

The distribution function / is related to the density 
per spin state by Eq. (|24|) and we assume that, as in 
equilibrium, the distribution functions for the two spin 
states are the same, i.e., f± = f-\ = f. This is true 
if the trap potential and the excitation operator of the 
collective mode are spin independent. 

The functional /[/] appearing on the rhs of Eq. (|3"Tj) 
is the collision integral. It describes collisions between 
atoms with opposite spin and depends on the differential 
scattering cross section as [33] 

I[f}= J d 3 Pl J d^v-v^fhil-Af'Xl-Afi) 

-/'/((l-A/Xl-A/r)], (34) 

where p and p 1 are the incoming momenta, p' and p[ 
are the outgoing ones; v and vi are the incoming veloc- 
ities p/m and Pi/m, respectively, f2 is the solid angle 
formed by the incoming relative momentum p — p 1 and 
the outgoing relative momentum p'— Pi of the two atoms, 
/ = /(r,p,i), h = /(r, Pl ,t), /' = f(r,p',t), etc. The 
factors of the type (1 — Af) are absent in the classical 
Boltzmann equation. They are a consequence of Fermi 
statistics and ensure that an atom cannot be scattered 
into a state which is already occupied. This Pauli block- 
ing effect can result in a strong reduction of the collision 
rate, especially at low temperatures. 

In order to study the collective modes of the trapped 
gas, we consider a small deviation Sf = f — f eq of the 
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distribution function from the equilibrium one. Usually 
Sf is strongly peaked at the Fermi surface, but it can 
conveniently be written as 

6f(v, p, t) = f eq (r, p)[l - Af eq (r, p)]$(r, p, t) , (35) 

with a smooth function $ [33|. Expanding the Boltz- 
mann equation (|31[) to linear order in the deviations 
from equilibrium, and considering that a change Sf of 
the distribution function results in a change SU of the 
self-energy, we obtain 

f eq (l - Af eq )($ + £ ■ v r $ - v. r (y T + E/e,) ■ v p $ 

+ /3— • V r <5f/) - -J[$] (36) 
m / 

with the linearized collision integral 

![$]= / d 3 P l /" dO^|v-Vx|/ e Je ? l 

x (i_^)(l-^0(* + *i-*'-*i)- (37) 

Now we have to specify SU , which appears in the lhs of 
Eq. (|36[) . Since we neglect any possible momentum de- 
pendence of U and SU, it is clear that SU can be written 
as 

SU(r,t) = J d 3 py(p)Sf(r,p,t), (38) 

where 7(p) is the functional derivative SU[f}/Sf, taken 
at f eq . In the low-temperature limit, Sf is so strongly 
peaked at the Fermi momentum pp that "f(p) may be 
replaced by a constant 70 = i{pf), as in Fermi liquid 
theory (34J. This results in SU = joSp, with 

Sp(r,t) = J d 3 pSf(r,p,t) = J d 3 pf eq (\-Af eq )<5>. 

(39) 

In addition, by choosing a particular form for Sf, namely 
fif = df eq /dfi, we can identify 70 with the derivative 
dU eq /dp eq taken at constant T, i.e., 

BU 

5U{v,t) = -^ *p(r,t). (40) 

dp eq p e ,(r),T 

In the present work, we shall assume that Eq. (|4"0")) is 
a reasonable approximation also at higher temperatures, 
although it cannot be rigorously justified in this case. 

Equation (|36p together with (|4"0"|) constitutes the start- 
ing point for our study of collective modes with in- 
medium effects. It is a generalization of the Boltzmann 
equation used in Refs. [3, [la, E3 m the case ^ = 9P 
(Hartree approximation). 

B. Trial function 

As mentioned in Sec. lIII Al the function $(r, p, t) char- 
acterizing the deviation from equilibrium can be sup- 
posed to be smooth in phase space. This allows us to 



make a simple ansatz for $ with a small number of coef- 
ficients rather than solve the linearized Boltzmann equa- 
tion (f3"6"|) exactly. 

For any collective mode of interest, the trial function 
<&(r, p, t) has to contain at least those terms which are 
necessary to generate the velocity field u(r, t) character- 
izing the mode [35| . The presence of a velocity field mod- 
ifies f eq into 

/(r,p,t) = / eg [r,p-mu(r,t)] (41) 

and leads to a deviation 

Sf~-0f eq (l-Af eq )p-u, (42) 

i.e., the trial function $ must at least contain a term 
proportional to p • u. When this term is inserted into the 
linearized Boltzmann equation (I36p . the operator p/m • 
V r — V r (Vr + U eq ) ■ V p on the lhs of Eq. f3"6")) generates 
new terms, as do the SU term on the lhs and the collision 
term I on the rhs. In general, the number of terms is 
infinite and the system cannot be closed. 

However, in the case of an ideal gas (U = SU = 1 = 0) 
in a harmonic potential Vr, and if u is at most linear in 
the coordinates, it is possible to solve Eq. (|3"6"|) with a fi- 
nite number of terms. For instance, a term proportional 
to xp x generates terms proportional to x 2 and p\, and no 
other terms are needed. In the opposite limit of an ex- 
tremely strong collision term, i.e., in the hydrodynamic 
regime, a linear velocity field solves exactly the hydrody- 
namic equations if the equation of state can be approx- 
imated by a polytropic one, which is in many cases an 
excellent approximation Q. We therefore assume that 
also in our case it will be a good approximation to in- 
clude into <& only those terms which appear in the ideal 
gas case (of course the coefficients will change). 

To be explicit, we will focus on the scissors mode 
(S), the radial quadrupole mode (Q), and the breath- 
ing modes (B). In order to check the consistency of our 
model, we will also consider the Kohn mode (center-of- 
mass or sloshing mode, K). The velocity fields and the 
corresponding trial functions for these modes are given 
in Table |U Note that in the case of the breathing modes, 
the axial and the radial modes cannot be treated sepa- 
rately because they are coupled (although the coupling 
may be weak in very elongated traps). 



C. Frequency and damping of collective modes 

By inserting each trial function <!> into the linearized 
Boltzmann equation (|36[) and taking moments of the 
equation, namely multiplying it by any of the terms con- 
tained in $ and then integrating over r and p, one obtains 
a set of homogeneous linear equations for the coefficients 
Cj. The condition that the coefficient matrix determi- 
nant is zero yields an equation for the frequencies of the 
collective mode. 
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TABLE I: Velocity fields and corresponding ansatz functions $ for the different modes under consideration. 



mode 



trap frequencies 



u(r,t) 



$(r,t)e' 



sloshing 
scissors 
radial quadrupolc 



radial 1 
axial J 



breathing 



(K) arbitrary 

(S) L0 X > LOy ~S> LU Z 

(Q) Wj = U)y — LOr 3> CiJ z 

(B) Wa; = LJy = L0 r ~> LU Z 



oc (1,0,0) 
oc (y, -a;,0) 
oc (as,-y,Q) 
oc (x,y,0) 
oc (0,0, z) 



CiX + C2Px 
ClXy + ClXPy + C3VP: 



d(x - y ) + c 2 (xpx 
Ci(a; 2 + y 2 ) + c 2 z 
+c 4 zp z + c 5 {p 2 



CAPxPy 

-ypy) + c3(pl -pI) 

2 + c 3 (xp x +yp y ) 

Py) + C&P 2 Z 



Let us start by the center-of-mass oscillation of the 
cloud, known as sloshing or Kohn mode {K). In experi- 
ments, this mode is used in order to determine the trap 
frequency with high precision [f| , since it is known to be 
an undamped oscillation with the frequency of the trap, 
independently of the interaction [l9l . |20T | . It is an impor- 
tant test of the consistency of our method to check that 
this property is preserved. 

Multiplying Eq. (|36|) (with $ = c\x+C2p x ) by x andp x 
and integrating over r and p, we obtain the following sys- 
tem of equations [note that the collision term (|37|) on the 
rhs of the Boltzmann equation (|36p does not contribute 
since I[x] = I[p x ] = 0)]: 



-IU) 



C\ — iV|C2 = , 



^jVf — "3") Cl — iiomN-[C2 



0. 



where 



iV T = 



rPeq 



(43) 
(44) 



(45) 



denotes the number of atoms per spin state and 

C= I <l- i r>l- i ,,Jl, i a-Af„ l )r^l (46) 



is a constant depending on the interaction. When calcu- 
lating the determinant, we obtain lj — lu x , independently 
of the interaction, as it should be. Of course, analogous 
results are obtained for the sloshing modes in the y and 
z direction. 

If we repeat the same steps as before for the case of the 
scissors, quadrupole, or breathing mode, an additional 
complication arises from the fact that now the collision 
term on the rhs of the linearized Boltzmann equation 
(|36p gives a non-vanishing contribution. More precisely, 
only the terms in $ which are quadratic in momentum 
contribute, since I[riTj] = Ifapj] = for i,j = x,y,z. 
Using the symmetry properties of the explicit expression 
for I[piPj], one can furthermore show that the r^r; and 
rkpi moments of I[piPj] (i,j,k,l = x,y,z) vanish and 
only moments involving two momenta, 



Hjkl 



d 3 rd 3 pI\piPj}p k pi 



(47) 



survive. Using I[p 2 ] — 0, one can show that these must 
be of the form 

lijkl = Is(SikSji + 6u6jk - g<5y<5fcz) , (48) 



where Is is the moment which is relevant for the scissors 
mode, i.e., 



Is = I xyxy = J d 3 rd 3 pl[p x py]p x p y . 



(49) 



Some more details on how Is is calculated are given in 
Appendix [Bl Now let us define the corresponding relax- 
ation time t as 



1 



J d 3 d 3 pl[p x Py}p x Py 



Ws 



f d 3 r d 3 p f eq (l - Af eq )plp2 m 2 



(50) 



where Ekm denotes the kinetic energy (cf. Appendix [All. 
It is essentially this parameter which governs the temper- 
ature dependence of the mode frequencies and damping 
rates. The definition (|50|) is identical with that intro- 
duced in Ref. [7[ . 

Using the definition (|50|) , we can write the equation for 
the scissors mode frequencies in the following form (see 
Appendix [C] for more details) : 



\ ih ) + {lo 2 - lo\ c+ ){lo 2 - w g jC _) = 0, (51) 



where u>s,h and ws,c± are the frequencies in the hydro- 
dynamic {lot — > 0) and collisionless {lot — > oo) limits, 
respectively. The hydrodynamic frequency is given by 



J S,h 



(52) 



and does not depend on the interaction. In the colli- 
sionless limit, there are two modes with different fre- 
quencies, corresponding to rotational [u oc {y, —x, z)] 
and irrotational [u oc {y, x, 0)] velocity fields. In a non- 
interacting gas, these two modes have the frequencies 
In the interacting case, they are 



US,cl± 



LO„ ± LO,, 



changed to 



J S,c 



i± = (^+^)(l- X /2) 



± ^4^(1 - X + X 2 /8) + {cut + <j^/4 , (53) 

where \ is the interaction dependent parameter defined 
in Eq. (|X6|) . 

For the radial quadrupole mode, the equation for the 
frequencies has the form 

lLu{LO 2 ~UJ 2 Q . cl )-l-{LU 2 -LJ 2 Q . hd )=0 7 (54) 
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the hydrodynamic frequency is, again, independent of the 
interaction and given by 

^Q.hd = 2w r > ( 55 ) 

while the frequency in the collisionless limit depends on 
the interaction: 

^ ci =4^ 2 (l-x/2). (56) 

In the case U — gp (Hartree approximation), x reduces 
to 3Ei nt /2E pot , where Ei nt and E pot denote the inter- 
action and potential energies (cf. Appendix [X| . In this 
case, our limiting frequencies (in the hydrodynamic and 
collisionless limits: lot — > and lot — > oo) agree with 
those of Ref. [H, 

In the case of the breathing mode, we obtain two 
frequencies lob±, corresponding to the axial and radial 
breathing modes. The low- lying mode (wb_) corre- 
sponds essentially to a motion in the z direction (ax- 
ial breathing mode), while the high- lying mode (lob+) 
corresponds to a motion in the radial direction (radial 
breathing mode). The equation for the frequencies has 
the form 

■ I 2 2 w ,2 2 \ 

tLu(lO -LJ Bcl+ )(iJ -W B)C ;_) 

-\(u 2 -loI m+ )(u 2 -LO 2 bm _)=0. (57) 

The expressions for the limiting frequencies LOB,hd± and 
ujb,cI± are given in Appendix iDl 

IV. RESULTS AND DISCUSSION 

In this section we will present our numerical results 
for the scissors, radial quadrupole, and radial breath- 
ing modes. We will discuss the frequencies and damping 
rates as functions of the temperature for finite values of 
1/kpa and for the unitary limit (1/kpa = 0). The fre- 
quencies lo and damping rates T are determined by the 
real and imaginary parts of the solutions of Eqs. (fSTj) . 
([53)1. and (f5"7) . respectively. 

So far, most experiments have been done on resonance, 
i.e., for 1/kpa = 0. Two exceptions are the study of 
the radial quadrupole mode over the whole crossover re- 
gion by Altmeyer et al. [3| and of the scissors mode 
at 1/kpa — —0.45 by Wright et al. [6j. In addition to 
these two experiments, we will compare our results to the 
scissors, quadrupole, and breathing mode experiments at 
unitarity described in Refs. @, 0]- 

A. Radial quadrupole mode at 1/kpa — —1.34 

In the first experiment on the radial quadrupole mode 
on the BCS side of the BEC-BCS crossover [J], the trap 
had frequencies Lo r — 2ir x 370 Hz and lo z = 2ir x 22 Hz 



and contained N = 400000 6 Li atoms. The highest mag- 
netic field used in this experiment, corresponding to the 
weakest interaction, resulted in 1/kpa = —1.34, which is 
the value we will consider here. The temperature is unfor- 
tunately not known, but we assume that it was between 
0.03Tf, the lowest value ever reported by the Innsbruck 
group [361 ] . and Q.lTp, the upper value given in Ref. 0. 

Results for frequency and damping as functions of tem- 
perature are shown in the first two panels of Fig. [4l while 
the third panel shows the damping rate vs. the frequency 
This latter representation was proposed in Ref. [7] in or- 
der to get rid of the temperature, which cannot easily be 
experimentally determined. The single data point shows 
unambiguously the necessity of the inclusion of the mean 
field since the measured frequency (~ 2. lav) lies clearly 
above the limiting value for a collisionless gas without 
mean field (2u> r ). 

The theoretical curves shown in Fig. [4] represent differ- 
ent levels of approximation in the calculation. In order 
to see the effect of the different improvements of the the- 
ory, we include them one after another. We start with a 
classical Fermi gas (dotted lines), using Boltzmann dis- 
tribution functions (f eq = e -/3(p 2 /2m+v T -p )) m the cal- 
culation of the density profile, without any mean field 
effects and with the free cross-section, Eq. (f26|) . without 
Pauli-blocking factors (1 — Afeq) in the collision term. 
Within this approximation, the system shows hydrody- 
namic behavior (lo — > y/2u r ) at low temperature and col- 
lisionless behavior (lo — > 2lo t ) at high temperature, with 
strong damping Y in the intermediate regime. In the rep- 
resentation of r vs. lo, this results in a curve similar to 
a semi-circle. The hydrodynamic behavior at low tem- 
perature is of course an artefact of neglecting the Pauli 
blocking in the collision term and it is in clear contradic- 
tion to the measured frequency. 

In order to cure this problem, we include the effect of 
Fermi statistics (short dashes), i.e., we use the Fermi dis- 
tribution function = l/( e /3(p 2 /2m+V T -Mo) + 1) m fo e 

calculation of the density profile, and the Pauli-blocking 
factors (1 — Af eq ) in the collision term. At this stage, we 
still use the free cross section and we do not include any 
mean field. Due to the Pauli blocking factors, the colli- 
sion rate goes now to zero at low temperature, and there- 
fore the system approaches the collisionless frequency 2lo t 
in both the low and high temperature limits. The high- 
est damping, and as a consequence the lowest frequency, 
is reached at a temperature of ~ 0.3TV • Since no mean 
field is included, lo and T depend only on a single param- 
eter, namely on t [cf. Eq. ([54]) ]. Therefore the results 
lie on the same curve in the lo — T plane as in the case 
of a classical gas (dotted curve), but this time only the 
small part of the curve corresponding to large values of 
r is covered. Although the frequency at low-temperature 
is now in better agreement with the data, it is still too 
low, since we have not yet included the mean field. 

The third step consists in switching on the mean field 
U (long dashes), i.e., the density profiles are now cal- 
culated With f eq = l/( e P(P 2 /2m+V T + U-no) + and the 
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FIG. 4: Frequency and damping of the quadrupole mode for 1/kpa = —1.34. The experimental result is taken from Ref. 
[j. The three panels display (from left to right) frequency as function of temperature, damping as function of temperature, 
and damping vs. frequency. The different lines represent different levels of sophistication of the calculation: Starting from a 
calculation for a classical gas without any mean field and with the scattering cross-section in vacuum (dotted lines), we include 
the Pauli principle in the equilibrium density profile and in the collision integral (short dashes), then on top of that the mean 
field U (long dashes), and finally also the in-medium scattering cross section (solid line). The dash-dotted line represents the 
hydrodynamic frequency, wq,m = y/2u! r - 



U eq and 5U terms are included in the Boltzmann equa- 
tion (j36j) . We still keep the free scattering cross section. 
Our calculation is limited to temperatures above 0.06Tp, 
corresponding to the critical temperature of the system. 
The mean field does not have a dramatic effect on the 
damping, but it increases the frequency, especially at low 
temperature. Frequency and damping at the lowest tem- 
perature are now both in excellent agreement with the 
measured values. Note that the inclusion of the mean 
field modifies qualitatively the curve T vs. u> shown in 
the third panel of Fig. 2J 

Finally, we replace the free scattering cross section by 
the in-medium one, Eq. (|27|) (solid lines). Unfortunately, 
the good agreement between the theoretical results and 
the measured frequency and damping at low temperature 
is deteriorated: The resulting damping is too high by a 
factor of two and the frequency gets shifted downwards, 
although not dramatically. However, as already men- 
tioned, the calculation is limited to temperatures above 
~ 0.06Tp (the critical temperature Tc), while it is pos- 
sible that the temperature in the experiment was lower 
(the presence of a small superfiuid region in the center 
of the trap would not contradict the observation of the 
collisionless frequency [II])- Extrapolating the damping 
curve obtained with the in-medium cross section to lower 
temperatures, it seems that the result obtained with in- 
medium cross section is not necessarily inconsistent with 
the experiment. Additional experimental data points at 
higher (and known) temperatures could help to settle this 
question. 



B. Scissors mode at l/kpa = —0.45 

Shortly after the quadrupole mode, the Innsbruck 
group studied the scissors mode at 1/kpa = —0.45 and 



at unitarity (1/kpa = 0) In this experiment, the 
trap had frequencies uj x — 2ir x 830 Hz, u y — 2tt x 415 
Hz, and ui z = 2tt x 22 Hz, and contained N = 400000 
Li atoms. The frequency u> and damping T were mea- 
sured for constant interaction strength as functions of 
the temperature. The experimental data for the case 
1/kFCL = —0.45, taken from Ref. Q, are shown in Fig. 
[5] together with various theoretical results. As in Fig. [4] 
we display, in addition to to and T as functions of the 
temperature, T as a function of u>. 

The meaning of the different curves is the same as in 
the case of the quadrupole mode discussed in Sec. IIV Al 
The most curious result, which has already been ob- 
served in Ref. [l8| for the case I/Ufcl = 0, is that already 
with classical statistics (dotted line) one can reproduce 
quite well the observed frequencies: At low temperature, 

the frequency is the hydrodynamic one ( \J ^ + — 

2ir x 928 Hz), and with increasing temperature it rises 
towards the collisionless frequency (u> x +o-> y = 2tt x 1245 
Hz). What is most surprising is that the agreement is 
very good at temperatures well below the degeneracy 
temperature Tf, where the classical approximation is not 
justified at all, while it fails at higher temperatures. In 
fact, the high-temperature behavior of the frequency is 
not reproduced by any of the theoretical calculations, 
which might be due to the importance of the anharmonic- 
ity of the trap at high temperatures. The agreement 
between the measured frequencies and those of a classi- 
cal gas is, however, purely accidental, as one can see by 
looking at the corresponding damping rates. The classi- 
cal statistics leads to a damping which is much too weak 
at low temperatures (dotted line). 

In fact, as in the case of the quadrupole mode dis- 
cussed in Sec. IIV Al the lack of Pauli blocking results 
in a high collision rate, leading to a perfectly hydrody- 
namic behavior. The inclusion of Pauli blocking (short 



11 




dashes) strongly reduces the collision rate and therefore 
increases the damping at low temperatures, resulting in 
a very good agreement with the measured damping rates 
(except near the peak at T/Tf ~ 0.15, which is proba- 
bly due to the superfluid- normal phase transition). Note 
that, since the interaction is much stronger now than in 
the case of the quadrupole mode discussed above, the 
collisionless regime is not reached, although at low tem- 
perature the frequency increases strongly towards the 
collisionless one. This increase of the frequency is not 
observed in the experiment, because at these tempera- 
tures the system is already in the superfluid phase and 
therefore its frequency stays close to the hydrodynamic 
one, even if the collision rate is low. This effect cannot 
be described in the framework of the simple Boltzmann 
equation which does not include superfluidity. But also 
at higher temperatures, the agreement of the frequen- 
cies obtained with Fermi statistics (short dashes) with 
the data is not as good as that obtained with Boltzmann 
statistics (dotted line). 

The inclusion of the mean field (long dashes) leads 
to a small reduction of the frequency, while the damp- 
ing is slightly enhanced, improving the agreement with 
the data. Since our calculation is limited to the normal 
phase, the curves are restricted to temperatures above 
~ 0.2Tp (the critical temperature Tq)- The frequencies 
are now well reproduced for temperatures above ~ 0.3Tp, 
while they are still slightly too high between ~ 0.2 and 
~ 0.3Tp. The damping is in excellent agreement with 
the data for all temperatures above Tq- 

Finally, the inclusion of the in-medium cross section 
(solid lines) leads to a big disappointment: The agree- 
ment with the data, in particular for the damping, is 
completely lost. The results are very close to those of 
the classical gas, similar to the findings of Ref. (2§| for 
the shear viscosity of the unitary gas and of Ref. [7J] for 
different collective modes at unitarity. The reason is that 
the enhancement of the cross section (cf. Fig. [2]) cancels 
the effect of Pauli blocking. 

Apparently the present theory has a fundamental prob- 
lem. Maybe the quasiparticle approximation made in 



Sec. Ill CI is too crude (although the QP density profile 
coincides very well with the NSR one): There might be 
important corrections due to energy and momentum de- 
pendence of the self-energy [32]. Even the validity of 
the Boltzmann equation itself might be questioned: The 
T matrix approximation can result in a long correlation 
time, leading to non-Markovian (memory) effects [37| . 

C. Collective modes in the unitary limit 

In this subsection we will finally show results for col- 
lective modes in the unitary limit (1/kpa = 0). We 
will again compare with experimental results obtained by 
the Innsbruck group, Ref. @ for the case of the scissors 
mode and Ref. [7j for the radial quadrupole and breath- 
ing modes. In the experiment on the scissors mode, the 
trap parameters were the same as those stated in the 
beginning of Sec. IIVBI Our theoretical results and the 
experimental data are shown in the first row of Fig. [51 
the representations are analogous to those in Figs. [Hand 

El 

The quadrupole mode was studied in an axially sym- 
metric trap having ui r = 2tt x 1100 Hz, ui z — 2ir x 26 
Hz and containing N = 600000 6 Li atoms '7}. The 
theoretical results and the experimental data are shown 
in the second row of Fig. [5J For the radial breath- 
ing mode, the trap frequencies were somewhat higher, 
namely uj r = 2n x 1800 Hz, lu z = 2n x 32 Hz, but the 
number of atoms was again N = 600000 [7]. The cor- 
responding theoretical results and experimental data are 
displayed in the third row of Fig. [5J 

In the case of the scissors mode (upper row of Fig. [5J , 
the qualitative features of frequencies and damping as 
functions of temperature are reasonably well reproduced 
by all approximations [even by the classical gas (dotted 
lines), as it happened already in Sec. lIVB"] . As in the case 
of the scissors mode at 1/kpa = —0.45 discussed in Sec. 
IIVBI the data stay near the hydrodynamic frequency at 
low temperature because of superfluidity, which is not in- 
cluded in our theory. For temperatures above Tq, which 
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FIG. 6: Same as Fig. [4] but for the scissors mode, the radial quadrupole mode, and the radial breathing mode (from top to 
bottom) in the unitary limit (1/fcpa = 0). The experimental data are taken from Ref. [(| for the scissors mode and from Ref. 
[?J for the radial quadrupole and breathing modes. The dash-dotted lines represent the respective hydrodynamic frequencies 
of each mode. 



within our theory is given by Tq ~ 0.3Tp, the frequencies 
are best reproduced by the calculation including mean 
field and in-medium cross section (solid line). Unfortu- 
nately, as it happened already in the case 1/kpa = —0.45, 
the strong enhancement of the in-medium cross section 
leads to a damping which is much smaller than the ex- 
perimentally observed one. The damping as a function of 
temperature is best reproduced by the calculation includ- 
ing the mean field but using only the free cross-section 
in the collision term (long dashes). However, if one looks 
at the plot r vs. u), both approximations result in curves 
which are close to the data. This means that both ap- 
proximations might be compatible with the data, if it 
turned out that the temperature axis was incorrect. 



In this context it should be mentioned that the tem- 
perature measurement in the experiment is not evident 
and in addition not completely independent of the the- 
oretical model used in the analysis. For instance, in the 
analysis of Ref. 0], the method introduced by Thomas 
et al. in Ref. [38| was used, which requires, among other 
things, the knowledge of a parameter (3 determining the 
"effective mass" in the unitary limit. As a consequence, 
in addition to the statistical error in the determination 
of the temperature, there could be a sizable systematic 
error coming from theoretical uncertainties. This obser- 
vation underlines the usefulness of the representation of 
r vs. u>, which is independent of the temperature. 

Let us now look at the radial quadrupole mode (sec- 
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ond row of Fig. [6]). In this case, we have to admit that 
the differences between the different approximations are 
much smaller than their deviation from the experimental 
data, i.e., all approximations fail to give a satisfactory 
description of the data. The mean field has only a very 
small effect (difference between short and long dashes), 
and the in-medium cross section (solid line) leads to an 
additional deterioration as compared to the free one (long 
dashes) . It seems as if in the experiment the continuous 
transition from the hydrodynamic to the collisionlcss case 
happened at a much lower temperature than in any of 
the theoretical results. But also in the u> — T plot, where 
the experimental data follow a very well defined curve, 
the theoretical results are quite far from the data. Fur- 
ther studies within an improved theoretical framework 
are needed. 

Concerning the radial breathing mode (third row of 
Fig. , it is difficult to draw firm conclusions from the 
figures. As it was the case for the scissors mode, the 
results obtained with mean field but without in-medium 
cross section (long dashes) give a satisfactory description 
of the damping between Tc ~ 0.3Tf and ~ 0.8Tf, while 
the damping obtained with the in-medium cross section 
(solid line) is too weak. Because of their strong scat- 
tering, the frequency data seem to be compatible with 
any of the theoretical curves, except for the sudden rise 
of the frequencies obtained with mean field (long dashes 
and solid line) if one aproaches Tc from above. This 
rise of the theoretical frequencies does not come from a 
reduced collision rate, leading to collisionless behavior, 
but from a sudden increase of the interaction-dependent 
parameter \' (defined in Appendix |D|) . This effect is an 
artefact of our QP approximation since it occurs when 
the density profile gets flat at some point different from 
r = 0, which happens only with the QP densities close 
to Tc, but not with the NSR densities. 



V. CONCLUSIONS 

In this paper, we studied the collective modes of a cold 
Fermi gas with attractive interaction in the framework 
of the Boltzmann equation with in-medium effects. Our 
starting point was the T matrix approximation (ladder 
resummation) , which is known to give a satisfactory de- 
scription of the BCS-BEC crossover. We discussed the 
corresponding single-particle self energy and in-medium 
scattering cross section above the critical temperature 
Tc- Within the QP approximation, which we have to 
make in view of the Boltzmann equation, the self-energy 
acts effectively like an attractive mean field, but it avoids 
the pathological problems of the Hartree approximation. 
Although we used at the present stage the simplest pos- 
sible version of the QP approximation (neglecting energy 
and momentum dependence of the self-energy), it repro- 
duces well the equilibrium density profiles obtained with 
the NSR formula, except near unitarity close to Tc- The 
in-medium scattering cross section, which is also obtained 



from the T matrix, is strongly enhanced as compared 
with the free one at low temperature. This is a precursor 
effect of the transition to the superfluid phase. 

We derived the frequencies and damping rates of 
different collective modes (sloshing, scissors, radial 
quadrupole, and breathing modes) with inclusion of the 
"mean field" and in-medium cross section. To that end, 
we linearized the Boltzmann equation around equilib- 
rium. A first important result is that the frequency of 
the sloshing mode is equal to the trap frequency even in 
the presence of the medium effects, in accordance with 
the Kohn theorem. We stress that this result can only be 
obtained if one uses in the linearized Boltzmann equation 
the same mean field as in the calculation of the equilib- 
rium density profile. [It therefore seems dangerous to 
calculate the sloshing-mode frequency by inserting an in- 
teracting density profile (e.g., a measured one), into an 
equation like Eq. (B2) of Ref. Q derived from the Boltz- 
mann equation of a non-interacting gas.] 

The frequencies and damping rates of the different col- 
lective modes were evaluated numerically as functions of 
the temperature and for different values of the scattering 
length (parameter 1/kpa). Because of the experimental 
uncertainty concerning the temperature, we discussed in 
addition the "hydrodynamic circles", i.e., the damping as 
a function of the frequency. It was shown that using clas- 
sical statistics one clearly cannot reproduce the observed 
collisionless behavior of the modes at small temperatures. 
By the inclusion of Fermi statistics (Pauli blocking), this 
problem is solved. Including the mean field on top of 
Fermi statistics, which is the main topic of this paper, 
has only a relatively small effect, but nevertheless it leads 
to a significant improvement of the measured frequencies 
and damping rates. For instance, we can for the first 
time quantitatively explain the observed upwards shift 
of the quadrupole mode in the collisionless normal phase 
at weak coupling (1/kpa = —1.34), which has attracted 
a lot of attention. However, in the strongly interacting 
cases, in particular in the unitary limit, the frequencies 
depend much more strongly on the collision rate than on 
the mean field and the mean field effects are of minor 
importance. 

On top of that, we included the in-medium effects in 
the cross section, which determines the collision rate (re- 
laxation time). Since the in- medium cross section be- 
comes very large when one approaches the critical tem- 
perature, this compensates the effect of Pauli blocking 
and reduces the result of the full calculation to some- 
thing comparable with the result obtained for a classical 
gas, as already noted in Refs. [7|, [28|. Clearly, this is an 
up to now unsolved problem and needs further examina- 
tion. A first possible extension of the present theory is to 
take into account the energy and momentum dependence 
of the self-energy in a more involved QP approximation 

& 

In addition, the method of taking moments of the 
Boltzmann equation, although it proved very successful 
in the past, does not correspond to a full solution of the 
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Boltzmann equation. In particular, this method is in- 
sufficient for a description of the damping due to the 
anharmonicity of the potential (trap and mean field) . In 
order to improve the calculation in this respect, a more 
involved numerical method for solving the Boltzmann 
equation is necessary. A very common method is that 
using so-called "test-particles" ; in the context of trapped 
fermionic atoms it was used, e.g., in Refs. fl5l. Il6j. 

Finally, a shortcoming of our present calculation is its 
limitation to temperatures above Tc, whereas the most 
interesting experimental results (e.g., extremely strong 
damping at a certain temperature) are probably related 
to the transition to the supcrfluid phase. At least in the 
weakly coupled regime, it should be possible to include 
the effects of superfluidity following the approach of Refs. 
PS EH- Work in these directions is in progress. 
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APPENDIX A: VIRIAL THEOREM 

Let us derive a relationship between equilibrium quan- 
tities which is very useful for reducing the number of 
interaction dependent parameters in the explicit expres- 
sions for the collective-mode frequencies. 

From the Boltzmann equation in equilibrium one ob- 
tains the following property of the equilibrium distribu- 
tion function: 

(^■V r -V r (V T + U eq )-V p )f eq (l-Af eq ) = 0. (Al) 

Notice that this property was already used in order to 
derive Eg . (1551) . 

In Ref. [17|, where the case U = gp is considered for 
the mean field, it is shown that multiplying equation ()A1[) 
by xp 2 p 2 and integrating over r and p, one obtains the 
virial theorem E kln - E pot + 3E mt /2 = S3), with 



E km = 2 / d 3 rd 3 p^-f eq , 
2m 



E pot = 2 / d 3 rV T p e 



E lnt =g I d 3 rp 



and 



eq ' 



(A2) 
(A3) 
(A4) 



Here we want to show that the same can be done for 
any function U eq (v) = U[p(r),T] in order to obtain a 
generalized virial theorem. 

As mentioned in the end of Sec. Ill Ft we define the 
rescaled coordinates fi = riUi/uj, in terms of which the 



trap potential reduces to Vr — muj 2 f 2 /2, and, conse- 
quently, also the density p eq and the mean field U eq be- 
come spherically symmetric in these coordinates. 

It is then found that the generalized virial theorem is 



E kin -E pot -Jd*r Peq 



~ dU eq 
or 



(A5) 



In the case U eq = gp eq (Hartree approximation), the last 
term can be integrated by parts, and the well-known re- 
sult for the virial theorem is recovered. 

In the general case, let us define the parameter \ char- 
acterizing the strength of the interaction as 



1 



X 



,dU e , 



E pot .J a 



Then the virial theorem can be written as 

Em 



Epot 



i-x- 



(A6) 



(A7) 



APPENDIX B: COMPUTATION OF THE 
RELAXATION TIME 

The equation determining the frequencies of the scis- 
sors mode contains the parameter r, defined in Eq. (|50p L 
Its evaluation is quite involved and we follow closely Ref. 
[ill ] in order to reduce the number of integrals. The in- 
tegral Is entering the definition of r can be most con- 
veniently computed if one observes that 



Is 



10 ^2 7 w 



(Bl) 



The explicit expression for Ig reads now 



1 f ^ ,3 ,3 , d<J IP- Pil 
— d rd pd p 1 dtt— 

10 / d\l m 



xf eq f eql (l-Af' eq )(l-Af' eql ) 
x [P 4 + (P- Pl ) 2 -(P-P') 2 



P-Pi) 2 ]- (B2) 



In order to reduce the number of integrals, one first in- 
troduces the variables k = p + p 1; q = (p — p 1 )/2, and 
q' = (p' — p)/2 (remember that |q| = |q'|). In terms of 
these variables, the factor in the second line of Eq. (|B2|) 
becomes 



1 



1 



AA 2 cosh P(E - po) + cosh /3k • q/2m 

1 



cosh P(E — po) + cosh /3k ■ q'/2m 



. (B3) 



with E = k 2 /4m + q 2 /2m + V T + U. The factor in the 
third line of Eq. §B^) reduces to 2q 4 - 2(q • q') 2 + (k • 
q) 2 /2 — (k • q') 2 /2. Note that the last two terms do not 
contribute to the integral since they are antisymmetric 
with respect to the interchange q <-> q'. Let us now 
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denote by 9,(f> and 9', ft the zenith and azimuth angles 
of q and q', respectively, with respect to k. The integrals 
over (f> and ft , which appear only in the third line of Eq. 
(|B2|) . can be done analytically and, writing 7 = cos 9 and 
7' = cos 9', we are finally left with 



1 f°° of 00 9 f°° 7' 
/ dff 2 \ dkk 2 / dqq 7 - 
t rn J Q J J < 



• da 

2Q-K 2 m Jo Jo Jo * ^ ^ 

x / dl I dl '( i+21 ' 2 ~ 3 ^V 2 ) 

1 

^ 

cosh /?(£? — fio) + cosh f3kq r y/2m 
1 

V 

cosh/3(£' — /x ) + coshf3kq r y' /2m 



(B4) 



This five-dimensional integral is evaluated numerically 
using a Monte-Carlo algorithm. 



APPENDIX C: FREQUENCIES OF THE 
SCISSORS MODE 

In order to determine the frequency of the collective 
modes, one has to take moments of the Boltzmann equa- 
tion. It is thus useful to rewrite Eq. (|36| ) as 



f eq (l - Af eq )$ - (V p f eq ) ■ V r T$ 



dU, 



eq 



dp 



a; 



Sp) 



+ T{V r f eq ) ■ V p $ = -/[$] , (CI) 

and to denote the three terms on the lhs by (i), (ii), 
and (iii). When integrating over r and p, the following 
identities arc useful: 



V p / eg = feg(l - Af eq )p, 

m 

dU R 



9fe9 U„^\ Af eq )h 1 



d_ 

dpi 



1 d 

mQ 2 f dr 



(C2) 
(C3) 
(C4) 



As an example, we report in more detail the derivation 
of the modified frequencies of the scissors mode. The 
corresponding trial function <3> is given in Table |TJ 

The contributions of the term (i) to the moments of 
Boltzmann equation are: 



d r d p (i) xy ■■ 



iuE pot (l + (fx) 



3f3m 2 L0 2 u! 2 



CX 



d 3 r d 3 p (i) xp y 
d 3 r d 3 p (i) yp x 

d 3 rd 3 p(i)p x p y = - 



- - ofl 9 13 ° 2 

3pw 2 

iuE pot m 

iujm 2 E kin _ ibJt 
3/3 6 C4 ' 



(C5) 
(C6) 
(C7) 
(C8) 



where the parameter ipx is defined as 



1 



i3 ~ ~2 3Pe q dU eq 

a rr 



1 



5E pot J dr dr i + l.™? 

mto z r or 

The density variation is given by 

X -rut 1 d Peq 1 

op = e exxy 



/3 3/i 

df-t 



(C9) 



(CIO) 



Using this, we can write for the contributions of (ii) to 
the moments of the Boltzmann equation: 



d 3 rd A p {ii)xy = 0, 



(Cll) 



e- iut cx , (C12) 



d 3 rd 3 p{ ll )xp y = {l+ f l)E 2 VOt ^ t 

d 3 r d 3 p (it) yp x = {1+ J l)E 2 pot e-^cx , (C13) 

d 3 rd 3 p^)p x p y = !^» e --'( C2 + C3 ) . (C14) 

Finally, the (iii) contributions to the moments of 
Boltzmann equation are: 

/ d 3 rd 3 p(iii)xy=-^e-^(^ + ^), (C15) 

J opm y 

d 3 rd 3 p(iii)xp y = - T ^^e- iut c 4 , (C16) 

d 3 rd 3 p (iii)yp x = -T^-e-^CA, (C17) 



3/3 

d 3 r d 3 p (iii) p x p v = . 



(C18) 



With these results, the equations for the coefficients Cj 
read 

iuj{\ + ifx)cx + mu}yC 2 + mWjC 3 = , 
(1 + <px)cx — irriLOC2 - (1 - x)m 2 w^C4 = , 



(1 + ipx)cx - imwcz - (1 - x)m 2 ujlc4 = , 



Co '+ (•;; -r iiU )<■] = I). 



(C19) 



The system has a solution if 

ito 



J - (col + co 2 y )} + 2w2 (^ + - X/2) 

+ (w*-w*) 2 (l-x) = 0, (C20) 

which is Eq. (fSTj) with ws,fui and u>s,d given by Eqs. ([5 
and 



APPENDIX D: FREQUENCIES OF THE 
BREATHING MODES 



We consider a trap with frequencies lo x = io y = io r and 
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2 | 2 

x +y 



xp x + yp y 



zp z 



2 , 2 
Px + Py 



2 i 2 

x- +y 



XPx + VPy 



ZPz 



Px+Py 



2(l + 2 yi -y 3 ) 



-y>3) 



m^u;^ (1 — x) 



2m' 2 u' 2 . 

3^(l + yi) 
5 — 5 

yi-¥>3 

m 

2+3^i — tps 
m 

itu 

2m^j!(l-x) 

2m^J^(l- x ) 



A 2 





-2m^(i_2 X + 2 X ') 
2m^ 2 ( X -2 X ') 



2^-37 



iw + 



m^( x _ 2x ') 
-2m^(l-| X + x') 



iuj I 1 

2 "r" 3t 



2 3r 



In Table IdI each line is obtained by taking one moment 
of the Boltzmann equation. For example, the third entry 
in the first column (1/m) is the coefficient in front of 
C3 if Eq. ()Clj) . with $ as given in the last line of Table 
HI is multiplied by x 2 + y 2 and integrated over r and p. 
Notice that the coefficient A appears only trough the zp z 
term. The new interaction dependent parameters x' an d 
<P3 that enter in the table are defined as: 



X 



2E, 



pot 



(f3 



1 f,3~~ dJJ eq 

a rrp eq 



i 



dr i-| 1 ^m. 

muj z r or 



(Dl) 
(D2) 



where \' reduces to 3Ei nt /2E po t in the Hartree case. 
From the determinant of the 6x6 matrix given in Table 
|D| we obtain the equation for the frequencies, Eq. (f57l) . 
The frequencies in the collisionlcss limit are given by 



The hydrodynamic frequencies are 



J B,hd± 



■ CO. 



c ± Vc 2 + d 



with 



r 3[16 + 25yi - 25x(l + Vi)] 



(D6) 



c =25 X 2 (2 + A 2 )(l + tpx) + 6x'(2 + A 2 )(8 + 5</>i) 
+ (5 + 4A 2 )(16 + 25<^i) - 24(2 + A)Vs 
+ x[-157 - 175<^i - 30x'(2 + A 2 )(l + tp{) + 48^ 3 
+ A 2 (-116- 125(^1 + 24^)], (D7) 

d = - 36A 2 [16 + 25(^1 - 25x(l + tpi)] 

x [32 + 50(^1 + 25x 2 (l + tpi) + 6x'(8 + 5ipi) 

- 3x(22 + 25</?i + 10x'(l + <pi) - 8cp 3 ) - 24<^ 3 ] , 

(D8) 



± Vo 2 



J B,cl± 



16 + 25^i - 25x(l + (pi) 



(D3) 



with 

a =25x 2 (l + A 2 )(l + Vl ) + 2[x'(2 + A 2 )(8 + 5<pi) 
+ (1 + A 2 )(16 + 25^i) - 4(2 + A 2 )v3 3 ] 
+ x[-10x'(2 + A 2 )(l + (px) - 3(1 + A 2 )(22 + 25<pi) 
+ 8(2 + A 2 )H, (D4) 

b = - 4A 2 (2 - x)[16 + 25^ - 25 X (1 + <pi)] 

x [32 + 5(Vi + 25x 2 (l + (pi) + 6x'(8 + 5(pi) 

- 3x(22 + 25^i + 10x'(l + (pi) - 8<ps) - 24^ 3 ] , 

(D5) 



Note that since b and d are proportional to A, the low- 
lying limiting frequencies (w| cl _ and w\ hd-)i corre- 
sponding to the axial breathing mode, tend to zero in 
the limit of a very elongated trap (A — > 0). 

In absence of the mean field, x = x' — (fi = (f3 = 
and the frequencies reduce to the known expressions 



J B,cl± 



2uj 2 r {\ + A 2 ± y/l - 2A 2 + A 4 ) , 



(D9) 



J B,hd± ., 



= -f (5 + 4A 2 ± ^25 + 16A 4 - 32A 2 ) . (D10) 
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